Thermal Conductivity Tensor in YBa2Cu307_2,: Effects of a Planar Magnetic Field 
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We have measured the thermal conductivity tensor of a twinned YBa2Cu307_i, single crystal as 
a function of angle 6 between the magnetic field applied parallel to the Cu02 planes and the heat 
current direction, at different magnetic fields and at T=13.8 K. Clear fourfold and twofold variations 
in the field-angle dependence of Kxx and n^y were respectively recorded in accordance with the d- 
wave pairing symmetry of the order parameter. The oscillation amplitude of the transverse thermal 
conductivity found to be larger than the longitudinal one k.^^ in the range of magnetic field 

studied here (0 T < _B < 9 T). From our data we obtain quantities that are free from non-electronic 
contributions and they allow us a comparison of the experimental results with current models for 
the quasiparticle transport in the mixed state. 
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PACS numbers: 74.25.Fy, 74.72. Bk,72.15.He 



I. INTRODUCTION 

The processes which influence the thermal transport of 
the quasiparticles (QP) under the application of a mag- 
netic field and below the superconducting critical tem- 
perature are a matter of current research in the physics 
of high-temperature superconductors (HTS). One of the 
difficulties to obtain clear evidence for one or other mech- 
anism from thermal conductivity measurements that in- 
volves the QP, is basically related to the separation of 
the QP contribution from the measured thermal con- 
ductivity. Due to the fact that the experiments mea- 
sure the total thermal conductivity given by the sum of 
the contributions due to phonons, kP*^, and to the QP, 
K°\ special methods are necessary to separate the rela- 
tively large phonon contribution to the thermal transport 
in HTS. Basically two methods for this separation have 
been treated in the literature. One method is based on 
a phenomenological description of the field dependence 
of the electronic contribution of the longitudinal thermal 
introduced first by Vinen et al.EJ 



conductivity k°}.^{T- 

and used in Refs. B'BH to estimate k°'j.(T). The other 
method, presented by Zeini et alJj, is based on simulta- 
neous measurements of the longitudinal and transverse 
thermal conductivity and the assumption of a field inde- 
pendent Hall relaxation time. 

To study thermal transport properties which depend 
on the QP contribution we have chosen to measure the 
thermal conductivity tensor as a function of angle 9 be- 
tween heat current direction and the magnetic field ap- 
plied parallel to the Cu02 planes at different magnetic 
fields. Because the oscillation amplitude of the thermal 
conductivities as a function of angle is related only to the 
QP contribution, most of our results are basically free 
from phonon contributions and can help to test scatter- 
ing models for the QP. We discuss our results below in 
terms of two mechanisms that influence the behavior of 
QP in the mixed state of HTS: Andreev scattering (AS) 
is currently discussed ii^n^he literature as possible scat- 
tering mechanism of QPcl'B and the Doppler shift (DS) as 
the main effect that changes the energy spectrum of the 



QP with fieldBBfil. 

The infiuence of the scattering mechanism on the ther- 
mal properties depends also on the symmetry of the order 
parameter. In the last years a large number of measure- 
ments has confirmed the d— wave pairing symmetry for 
most of the HTS. Due to the d— wave symmetry of the su- 
perconducting order parameter in HTS, extended quasi- 
particle states in the mixed state - associated with a gap- 
less structure at the nodes - are expected to contribute to 
the thermal transport. These measurements have been 
described in terms of a DS in the energy spectrum of QP 
and the scattering of QP by vortices (AS) taking into 
account the presence of line nodes on the Fa'mj surface, 
i.e., within a pure (i3;2_j^2 -pairing symmetry! Oli3. 

Try-crystal phase-sensitive measurementsEj have de- 
termined that the pairing symmetry of the order pa- 
rameter is dx2_y2, which predominates with a small, if 
any, imaginary component (< 5%) id^y at T < Tc- Al- 
though symmetry experiments should not be significantly 
affected by the development of a small idxy contribu- 
tion, the quasiparticle transport might show (specially at 
T <C Tc) quantitative differences from the single d^^-y^- 
pairing state due to the absence of line nodes in the 
resulting superconducting dx2_y2 + i^Ly-gap. This new 
state proposed originally by Laughliril3 could arise either 
by decreasing teujipeMture or increasing magnetic field as 
argued in Refs.oEJlla. Up to now, however, there is not 
enough experimental evidence for such a state although 
results on the magnetic field dependepee-pf the longitu- 
dinal thermal conductivity in BSCjCpliZHla and tempera- 
ture dependence of the Hall angletj could be explained 
taking into account the violation of both parity and time- 
reversal symmetry. In this article we recover this idea 
and attempt to establish the highest value, if any, for a 
component id^y of the superconducting order parameter 
within the current models for the thermal transport in 
the HTS cuprates. Our results agree with the dx2-y2- 
pairing symmetry of the order parameter and provide a 
maximum value of ~ 35% for the component idxy 

From the phenomenological side our work provides in- 
teresting details. Namely, we show explicitly how the 
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transverse thermal conductivity Kxy can be calculated 
from the experimental transverse temperature gradient 
VyT and the longitudinal thermal conductivity Kxx for 
magnetic fields applied parallel to the Cu02 planes. As in 
the case of a perpendicular magnetic field, the resulting 
Kxy depends on the ratio of the transverse temperature 
gradient to the longitudinal one and represents a physical 
properUijOf the crystal involving only contributions from 
the QPE3. We also show that the oscillation amplitude of 
the transverse thermal conductivity, n'^y, is larger than 
the longitudinal one, kJ^^., in the whole magnetic field 
range studied (0 T < B < 9 T). From the extracted elec- 
tronic contribution of Kxx we attempt also to calculate a 
Hall-like angle and show that it increases with magnetic 
field in qualitatively agreement with the AS model. 

Our article is organized as follows. In the next section 
we describe some experimental and sample details; in 
Sec. [I] we describe and discuss our results and is divided 
in four subsections. A summary is presented in Sec. |^ 



II. EXPERIMENTAL AND SAMPLE DETAILS 

For all the measurements presented here we used a 
twinned YBa2Cu307_;r single crystal with dimensions 
(length (Z)x width (w)x thickness (d)) 0.83 x 0.6 x 
0.045 mm'|-ABd critical temperature Tc — 93.4 K, studied 
previouslyt3't2l. With polarized light microscopy we have 
determined the positions of the twinning planes. Ac- 
cordingly, the heat current direction was adjusted along 
the a/h axes which were parallel to the width or length 
of the sample. The use of a twinned crystal rules out 
the influence of orthorhombicity on the thermal trans- 
port properties studied in this work. In particular we 
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can assume that Kxx{T,B) 



^yy 



iT,B). We note that 



the influence of the orthorhombicity can be observed in 
untwinned crystals only if the inipiirity scattering rate 
is small enough. As shown in Ref.E2l, untwinned crystals 
with relatively large impurity scattering show similar an- 
gle dependences as twinned ones. 

The longitudinal (VxT) and transverse (VyT) tem- 
perature gradients have been measured at T = 13. SA" 
with previouslA2-|Calibrated chromel- const ant an (type E) 
thermocoupleiil and a dc picovoltmeter. Our definition 
of sample axes, temperature gcadients and field direction 
is similar to that used in Rcf.t3. If the heat current Q is 
along -I- J? and a positive 90° angle 9 on the (x, ?/)-plane 
is along -\-y, then there is a positive transverse thermal 
gradient for the magnetic field B at -1-45°, see Fig. |l|. 

The magnetic field was applied perpendicular to the 
c-axis. Special efforts were made in order to minimize 
the misalignment of the plane of rotation with the Cu02 
planes. We estimate the misalignment angle to be smaller 
than a ~ 0.5°. An in-situ rotation system enabled us the 
measurement of the thermal conductivity as a function of 
the angle defined between the applied field and the heat 
flow direction along -I- a;, see Fig. |^. We used an initial 
temperature gradient VxT ~ 241 K/m (A^^T ~ 200 mK 




FIG. 1: (a) Experimental arrangement, sample axes and ther- 
mal gradient definitions. The " -I- " and "-" signs of the ther- 
mocouples refer to the hot and cold ends when the magnetic 
field is applied at 45°. The longitudinal gradient along x is 
always negative. The misalignment angle a was minimized 
step-by-step measuring the angle dependence of the longitu- 
dinal thermal gradient till reaching a satisfactory symmetrical 
curve, (b) Top view of the sample arrangement with a sketch 
of the d— wave symmetry of the superconducting order pa- 
rameter. 



for our crystal) &i 6 = 90°, i.e. with the magnetic field 
perpendicular to the heat flow, and we recorded angle 
scans and magnetic field dependeiit;e-.with a field-cooled 
procedure to avoid pinning effectal3'Ea. That means that 
the angle and magnetic field scans at constant tempera- 
ture were done in such a way that after taking a reading 
the sample was driven to its normal state by heating up to 
a few Kelvin above the superconducting critical temper- 
ature and immediately was cooled down again to 13.8 K 
at constant field. We will show in Sec. Ill C experimental 
evidence that demonstrates how pinning of vortices influ- 
ences the field dependence of the measured oscillations, 
even for the case of parallel field used in this work. 
The possible infiuence of demagnetization and fiux pin- 
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FIG. 2: (a) Field-angle dependence of the normalized tem- 
perature gradient (A^r — 200)/Z at 13.8 K for a magnetic 
field strength of 8 T (solid circles) and 3 T (open circles). 
The heat flow was constant during the whole scan and chosen 
to produce 200 mK of temperature difference at the initial- 
scan angle 6 = 90°. Dashed lines are fits to the function 
a(cos(4^) — 1) with a(T,B) as free parameter. 

ning effects associated with the rectangular shape of the 
sample, see the commept, and reply in Ref.a, has been 
checked by Aubin et al.E3. The authors concluded that 
the relative orientation of vortices and the heat current, 
and not the sample geometry, governs the angle depen- 
dence. In our case, due to the relative large applied fields 
in comparison with the lower critical field Hd and the 
field-cooled procedure, an influence of the magnetization 
and flux pinning on the angle dependences can be ex- 
cluded. 



III. RESULTS AND DISCUSSION 

A. Angle dependence of the longitudinal n^x and 
transverse K^y thermal conductivities 

Figure || shows the field-angle dependence of the tem- 
perature gradient (A^X-200) /I at 13.8 K for 8 T and 3T. 
As pointed out in Ref .E3 the longitudinal thermal conduc- 
tivity Kxx oc Q/A^T and hence, the longitudinal temper- 
ature gradient, shows fourfold symmetry with minima 
(maxima for the thermal conductivity) along the nodal 
directions of the superconducting gap, see Fig. ||(a). We 
note that this effect represents a strong evidence for the 
d2,2_j^2-pairing symmetry of the order parameter, con- 
firming the presence of a large number of quasiparticles 
at those particular orientations. 

The curves in Fig. || can be fitted satisfactorily with the 
simple function a (cos (40) — 1) with only a as fit param- 
eter which depends on temperature and magnetic field. 
The fact that minima are found in the field-angle profile 
of the temperature gradient {A^T — 200)// at the nodal 



directions could be explained in terms of Andreev scat- 
tering of QP by vortices. According to this effect if at 
any point the quasiparticle energy, as viewed from the 
superfiuid frame, equals the gap then, at that point, An- 
dreev refiection occurs: the quasiparticle is converted to 
a quasihole reversiiia; its velocity and its contribution to 
the heat transports. This mechanism of thermal resis- 
tance is always induced when spatial variations of the 
amplitude or phase of the order parameter take place. 
Since the QP in the mixed state are in the presence of a 
phase gradient (superfiuid flow) they should also have a 
DS in their energy spectrum given by the scalar product 
of the momentum and superfluid velocity, p • Vg as viewed 
from the laboratory frame. That phase gradient is then 
the necessary mechanism that induces the Andreev re- 
flection. Thus, when a quasiparticle with momentum p 
is moving parallel to the magnetic held no DS occurs and 
hence, no Andreev reflection. In other words, the field 
acts as a filtero for that quasiparticle that contributes to 
reduce the total temperature gradient. For different ori- 
entations of the quasiparticle momentum with the field, 
there is a finite probability for the Andreev reflection 
and thereby, for an enhancement of the thermal resis- 
tance. The directionality of this scattering mechanism 
is therefore the responsible feature for the identification 
of the preferential directions of the QP by rotating the 
magnetic field parallel to the basal planes. 

On the other hand, a DS in the energy spectrum of the 
quasiparticle could also give rise to fourfold symmetry 
in the longitudinal thermal conductivity by rotating the 
field parallel to the Cu02 planes, but in principle, with 
opposite sign. The usual effect would be to produce an 
excess of quasiparticles in the direction perpendicular to 
the field, and |tii|eceby, reducing the thermal resistance in 
that directionaEa. However, the DS affects both the car- 
rier density as well as the scattering rate of QP, and since 
the latter dominates at high enough temperatures the QP 
move easier parallelptp the field and the sign turns out 
to the correct oneEjCJ. We note that, although it is not 
clear which effect (AS or DS) should prevail at low fields, 
the quasiparticle mean free path shortens with increas- 
ing field as might be expected if the vortex are moving 
closer together and therefore, AS should dominate the 
heat transport at high enough fields. 

In Fig. 1^ we show the field-angle dependence of the 
transverse temperature gradient AyT/w at two magnetic 
fields. We note that the sign of this effect is just that ex- 
pected within the pure AS picture, i.e., the QP with mo- 
mentum parallel to the magnetic field does not suffer An- 
dreev reflection, and therefore, enhances the transverse 
temperature gradient in that direction. In other words, 
the temperature of a side of the sample is enhanced when 
the held is oriented to that side. As argued in Ref.Ej 
and since the transverse signal is expected to be free of 
non-electronic contributions in contrast to the longitudi- 
nal temperature gradient, this fact could be deflnitive in 
distinguishing between Andreev reflection and DS. How- 
ever, the inclusion of an impurity scattering rate in the 
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FIG. 3: Angle dependence at two magnetic fields, 8 T (solid 
circles) and 3 T (open circles), applied parallel to the Cu02 
planes of the (a) longitudinal and (b) transverse thermal con- 
ductivities. The continuous curves are obtained from Eq. (^) 
with the parameters given in the text and a pure d^2_y2—wave 
order parameter. 



E 

<l 




ign of both 
3. Further- 



calculations of the DS effect accounts f| 
fourfold and twofold oscillations as wel 
more, in the case of BSCCO it has been argued that at 
the range of magnetic field and temperature studied here, 
both effects are taking place, and depending on the val- 
ues of the parameters governing the gap at the nodes, 
i.e, vf/v\\ and kp, along with an impurity scattering rate 
Fq, different regimes in which DS dominates oyer An- 
dreev scattering or vice versa, could be IdentifiedEZI. Un- 
fortunately, there is no yet theory for both longitudinal 
and transverse thermal conductivities which takes into 
account both effects when the field is oriented parallel to 
the Cu02 planes. Also, the theory for the thermal trans- 
port for our case (i.e. parallel field, low temperatures 
and finite impurity scattering) under the influence of DS 
is still not yet developed. Therefore, in most of this pa- 
per we restrict ourselves to a quantitative comparison of 
our experimental results to the predictions obtained from 
the AS mechanism; deviations from the predictions may 
well indicate that the other mechanism could influence 
the thermal transport, as suggested in RefE3. 

For a better comparison with theoretical predictions 
we need to obtain the transverse thermal conductivity 
from the measured transverse temperature gradient. As 
in the case the field is applied normal to the basal planes, 



we calculate the transverse thermal conductivity, K^y as 
a function of the field angle 6 from the Onsager relations, 
using the fact that the net heat current flowing in the y- 
direction is zero. At constant temperature, this can be 
written as follows: 



\/xT{B,e) 
J\/yT{B,e 
" (VxT(i?))2 



(1) 



-45 45 
e (degrees) 



135 



where we have used the experimental fact that the total 
longitudinal temperature gradient (~ 241 K/m) is much 
larger than its fourfold oscillation (^ 3 K/m) as shown 
in Fig. 1^. This allows us to approximate {VxT{B, 9))^ ~ 
{yxT{B)Y in (j^) and thereby, the fourfold symmetry 
does not enter when calculating n^y. In other words, the 
transverse thermal conductivity has the same symmetry 
as the transverse temperature gradient but opposite sign. 
We note that (|l|) ensures us that nxy measures in fact a 
physical properties of the crystal as in the usual Righi- 
Leduc effect. Figure ^(b) shows the angle dependence 
of the transverse thermal conductivity at two different 
magnetic fields. 



FIG. 4; Field-angle dependence of the transverse temperature 
gradient AyT/w at 13.8 K for a magnetic field strength of 8 
T (solid circles) and 3 T (open circles). As in Fig. ^, the 
longitudinal temperature difference at the beginning of the 
scan (at 8 = 90°) was 200 mK. Dashed lines are fits using 
bsm{26) with b(T,B) as free parameter. 



B. Separation of the phonon and QP contributions 
'a la Vinen' and comparison with the 
Andreev-scattering mechanism 

In order to separate the phonon and quasiparticle con- 
tributions from the longitudinal thermal conductivity, 
needed for a quantitative comparison with the AS model, 
we have measured the field dependence up to 9 T with 
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FIG. 5: Magnetic field dependence of tlie longitudinal thermal 
conductivity Hxx with the field oriented along the direction of 
the heat current (0 = 0°) and along a nodal direction {6 — 
45°). The two continuous curves are obtained with (2) with 
n = 1(0.82) for the angle 6 = 0°(45°). 



the field oriented along both antinodes and nodes direc- 
tions. As shown in Fig. |^ the field profile can be rather 
well fitted with the empirical expressionEJ 



AT) 



(2) 



where, in general, j3e{T) is proportional to the zero- field 
electronic mean-free-path of the QP and n cajjuhe related 
to the nature of the quasiparticle scatteringtlBa. 

It is clear that the phononic and electronic contribu- 
tions at zero field, and k^'^, are independent of the an- 
gle 9. Further, we assume that Pe{T) remains unchanged 
by rotating the field parallel to the Cu02 planes at con- 
stant temperature. In this case the exponent n might 
be field-angle dependent due to the different types of ef- 
fective scattering which the QP experience at different 
orientations of the magnetic field. This assumption re- 
mains thereby valid either_when the main effect is the 
DS or the AS. As in Ref.cl we used the value n = 1, 
when the field was aligned parallel to the direction of the 
measured temperature gradient (parallel to the heat cur- 
rent, 9 = 0°) and it was allowed to take other values for 



different orientations of the field. As shown in Fig. ||, 
we get satisfactory fits to the experimental points with 
= 5.29 W/Km, k^J^ = 0.89W/Km, and n ~ 0.82, the 
latter when the field is oriented along the nodal direction 
6 = 45°. To the best of our knowledge, it appears that 
this is the first time such a dependence of the exponent 
n in (H) with angle 9 is reported. 

The microscopic model for the scattering rate of the 
QP by vortices using pAndreev refiection ri,(-B, p) was for- 
mulated by Yu et aln. The thermal conductivity can be 
calculated usinft|a-2D version of the Bardeen-Richayzen- 
Tewordt modelclEa, 



.el 



1 



' ^2p}'9aVgf3El ^^^^^ 



r(B,p,T) 



sech 



2fcsT 



(3) 

where r(B, p, T) may be taken as a relaxation rate given 
by the sum of the following scattering mechanisms act- 
ing in series: scattering of QP by impurities ri,iip(p), by 
phonons rph(p,T), by QP rpp(B,p,T) and AS by vor- 
tex supercurrents ri,(B, p, Tp. The parameter Ep is the 
QP energy, a and (3 denote the x or y directions on the 
plane of the sample, Vga is the group velocity along t 
a direction and c is the c-axis lattice parameter. In Ref 
an approximate expression for the effective rate of AS by 
vortex current was obtained which is basically given by 



r„ (X a„ ^ exp(-a^/(£;p), cj)) / ln{ay / ^)) , 



(4) 



In this expression the effective fiux-line-lattice parameter 
is given by oc Bc2/ B^, where Bc2 is the upper critical 
field parallel to the planes and an effective mass factor 
7 ^ (m* jra)^. The function / depends on the energy 
of the QP, and consequently on the gap A(p), and the 
angle between the direction of the momentum of t 
QP and the magnetic field. For more details see Ref. 
As we mentioned above, we would like to know how large, 
if any, can be the contribution of a term idxy in the order 
parameter compatible with our transport measurements. 
Therefore we use the following gap 



A(p,r) = A(r)(A, 



(5) 



which can be explicitly written taking into account the 
momentum components as 



cos(_H£)-cos(££) . 2sin(^)sin( a^) 
l-cos(2f^) l-cos(£f 



A(p,r) = A(T) ( -\ ^ ' +z i^.:;;;af ' i ■ (6) 



The percentage values indicated below for the id-j^y part means its weighted contribution relative to the total gap. 

Our measurements are done at constant temperature, scattering rates Pimp, Pph and Pqp do not depend on the 
therefore if we make the simplifying assumption that the momentum of the excitations, we can write Pimp + Pph + 



6 



Fqp = To ~ cte. At higher temperatures it would be 
necessary, however, to include a model for the energy 
dependence of the-iajslastic scattering rate Fqp as argued 
by many authorscElEJ. Nowadays there is consent that the 
peak observed in the longitudinal thermal conductivity as 
a function of the temperature is due to the competition 
between a decreasing inelastic scattering rate and the 
decrease in the number of quasiparticles with decreasing 
temperature. 

From (|) we obtain k^'^Fq = 7.2 x lO^^CW/Kms) at 
zero field assuming (i2.2_y2-pairing symmetry and maxi- 
mum gap A(r = 0) = Aq = 20 meV. Then, if we use 
the quasiparticle contribution obtained from the fits 
using (||), we get Fq^ ~ 11.12 ps. Using this value and 
the parameters from Ref.Q, i.e. Bc2 — 500 T, Ginzburg- 
Landau parameter k = 100, 7 = 4 and T~^{B = Bc2) — 
3/2F^^ (chosen to reproduce the experimental amplitude 
of the oscillation at 8 T) it is not possible to fit satisfac- 
torily the thermal conductivity HxxiB, 9) of Fig. ||(a) for 
both curves at 3T and 8T simultaneously. However, if 
we increase strongly the upper critical field to an unre- 
alistically high Bc2 — 9, 500 T or increase by a factor 
~ 5, Eq. (||) accounts for the fourfold oscillation observed 
in Kxx{B) (and hence, for the normalized temperature 
gradient {AxT — 200)//)) when the field is rotated from 
9 = +90° to 9 = -90° as shown in Fig. |(a). 

With the same parameters mentioned above, Kxy can 
be calculated from (||) as shown in the Fig. |^(b). We 
find agreement between the experiment and theory if, as 
for Kxx, we increase the Jpiervortex spacing 5 times 
the value defined in Ref.Qild. This result indicates that 
the AS rate needed to fit the experimental data is much 
smaller than the one obtained from the proposed raicro- 
scopic model and with the parameters used in Ref.cl. 

We can also extend the AS model including in (||) 
an additional component idxy to the pure dx2-y2-gap 
and compare it to the experimental angle dependence 
of Kxx{9), for example. The results shown in Fig. 6 indi- 
cate that the fourfold oscillation is no longer observable 
with a component larger than ~ 35% at 13.8 K. Because 
we compare the symmetry of both theoretical and exper- 
imental curves taking the normalized Kxx (see Fig. 6), 
this upper limit is not affected using other parameters in 
Eq. (|. 



C. Field dependence of the oscillation amplitudes 
kIj. and K^ly-. Doppler shift and the effect of pinning 
of vortices 
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FIG. 6: Angle dependence of the normalized longitudinal con- 
ductivity at a field of 8 T Kxx{6) / i^xx{0)- The data points are 
taken from Fig. ^(a) normalized by an arbitrary factor in or- 
der to compare with the theoretical predictions given by the 
continuous, dashed and dotted lines. These lines were calcu- 
lated with (^ assuming a pure d^2_y2 symmetry, plus a 35% 
and 50% idxy contribution added to the gap (see Eq. (^), re- 
spectively, and with the same parameters used to fit the data 
of Fig. I 



E~ 0.08 
^ 0.06 




FIG. 7; Magnetic field dependence of Kxx (open circles) and 
(solid circles). Solid lines are fits using Eq. (31). Dashed 
lines are fits using the DS model following Eqs. (3). The in- 
set shows the amplitude Kxy recorded with a non-field cooled 
procedure (*) and with a field-cooled procedure (solid circles, 
the same data as in the main figure). 



Other way to compare theory and experiment that 
does not need any empirical relation to separate the QP 
contribution is based on the measurement of the field 
strength dependence of both fourfold and twofold oscilla- 
tion amplitudes, which are in fact pure electronic effects. 
It is convenient to define the quantities 

K'iy = Kxy{9^0°)-Kxy{9^45°), 



= Kxxi9 ^ ~ Kxxi9 ^ 0°) , (7) 

which can be directly compared with any model that ac- 
counts for the contribution of the QP. In Fig. ^ we show 
the magnetic field dependence of n'^y and k^^. As shown, 
we find that the resolution of k^^. is smaller than that of 
K^y. The fourfold oscillation observed in Fig. |2| has a 
thermal background of ^ 200 mK, which makes more 



7 



difficult the observation of any small superposed signal. 
The thermal background in the transverse thermocou- 
ple due to a small misalignment in the x-direction was 
~ 8.5mK. 

We note that the oscillation amplitude of the trans- 
verse thermal conductivity K'^y is larger than the one of 
the longitudinal thermal conductivity k^^. Note the fact, 
that Eq. (||) with the parameters described above (includ- 
ing an intervortex spacing 5 times larger) accounts for 
this behavior, see Fig. |^. One would tend to believe that 
changing the effective mass factor 7 we may compensate 
the too high i3c2-value needed to fit the experimental 
data in Fig. |^. Note that 7 enters not only in but in 
other parts of Eqs. (||) and (^. Our numerical results 
show no agreement with the experimental data of Fig. M 



when we choose 7 > 4 changing Bc2 between 10,500 T 
and 500 T. 

With the proposed DS mechanism by Won and Makii^l 
that would inffuence the quasiparticle transport we can 
also get satisfactory fits to our experimental points either 
for n'^y or but not for both components with exactly 
the same parameters. As we noted above and in Ref.Ey, 
the predictions of the DS model for the thermal conduc- 
tivities with a parallel field were obtained in the super- 
clean limit (this means that no impurity scattering rate 
appears explicitly in the thermal conductivity tensor) 
and for the temperature range ksT » e = {l/2)\/vv'eB 
(v and v' are the Fermi velocities perpepdicular and par- 
allel to the Fermi surface, respectively )c3: 



2A 

ln(^^)ln(- 



4A 



(27r)2 A2 '1.75T 



/eH 



) sin 29 



77r2 T 2.2, , 2A 



1 . 2A , 4A , 

27r^ 1.75J A^ wWeH 
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(1 - cos (46*))] 



(8) 



where A is the energy gap and, in units of /ifcs, k„ = 
Tn/3Am, n the density of QP and m their mass. Both 
expressions provide the correct sign for the tensor com^ 
ponents. From these equations and as shown in Ref.c2l 
we obtain for the field dependence of the oscillation am- 
plitudes K^y = aBhi(l/bB) with a and b fit parameters, 
and k|^, = aB. These dependences also agree well with 
the experimental data, see Fig. 0. The fit of n'^y gives 
a ~ 0.015 (W/TKm) and b = 0.048T-1; from the ex- 
pression proposed for we get a ~ 0.010 (W/TKm). 
This discrepancy in a is not too large taking into account 
that neither the superclean limit used in developing this 
model nor the field-temperature range where the Eqs. (|^) 
are valid, refiect the situation studied here. Therefore, 
a more accurate treatment of the impurity scattering is 
necessary to check if the DS picture can fit quantita- 
tively the experimental data. It is striking that the DS 
description appears also to provide the experimental field 
dependences of the oscillation amplitudes. 

In the inset, we show the effect of the pinning of vor- 
tices on the oscillation amplitude k^^,. While the points 
taken with a field-cooled procedure show a monotonic 
behavior in the range of fields studied here, the points 
taken with a non-field-cooled procedure, i.e. in this case 
the field was increased from zero at constant temperature 
and always at the same initial angle, show a saturation 
above 3 T. This effect could have a different origin 
from the DS, and would indicate that vortex scattering 
becomes more important with increasing field. We note 
that the pinning of the Josephson-like vortices parallel to 



the planes is strongly affected by vortices perpendicular 
to the planes due to the misalignment of the crystal re- 
spect to the applied field. The results in the inset of Fig. |^ 
clearly stress the necessity of field-cooled experiments to 
get the true magnetic field dependence of thermal trans- 
port parameters. 



D. Hall-like angle and quasiparticles-dependent 
ratios 

The angle dependent conductivities we have measured 
allow us the definition of quantities which just involve 
contributions from the QP as in the usual Righi-Leduc 
effect. As noted above we do not have to assume strictly 
any method to separate electronic and phononic contri- 
butions from the total thermal conductivity to obtain 
quantities which can be directly compared with the the- 
ory of QP transport. As shown in the last section, some 
of these quantities are themselves the field dependence of 
both and However, the ratio n'^y/'^xx might pro- 
vide also an interesting result for comparing with theory. 
In fact, it describes an intrinsic property of the crystal 
studied, free of the experimental uncertainty to deter- 
mine the total power given to the sample. 

As shown in Fig. two regimes in the oscillation ratio 
k^^/k", can be distinguished: From zero to ~ 2 T the 
ratio decreases stronger with field than at fields above 
'-^ 2 T. An abrupt decrease of this ratio is predicted by 
Eq. (0) but it fits well the experimental data only above 
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FIG. 8: Magnetic field dependence of the ratio K%y/K%^ and 
f^xx/f^xy (inset). The continuous line is obtained with Eq. (^) 
and with the same parameters of Fig. ^ described in the text, 
whereas the dashed line is obtained with an upper critical 
field Bc2 = 650 T. 



~ 2 T, see Fig. g. This discrepancy at low fields could 
indicate that at i? < 2 T other mechanism influences the 
QP transport and partially overwhelms the influence of 
the proposed AS. In the same figure and to demonstrate 
the infiuence of the fiux-line-lattice parameter on the 
theoretical predictions we show the ratio (dashed line) 
calculated with the same parameters as before but with 
Sc2 = 650 T. 

In order to provide further analysis of our data, it is 
interesting to define a Hall-like angle as in the case of 
a normal field applied to the Cu02 planes, and to ex- 
plore its magnetic field dependence. Therefore, since in 
this case the angle 9 is an intensive parameter of the mea- 
surement along with the temperature T and the magnetic 
field strength B, it seems reasonable to perform an aver- 
age on the angle dependence of both conductivities over 
the range [—90, 90°] of 6. Also, we do not need to take 
into account the sign of K^y which indicates which side 
of the sample becomes hotter when the field is applied, 
characterizing the defiection direction of the QP. At con- 
stant temperature, we define the effective Hall-like angle 
as 



tan ail 



^ {\^^.yiB)\)e 



(9) 



where ( )g means the integration over 6 in the range men- 
tioned above. 

In Fig. ^ we show the field dependence of tanan in 
which we used the phononic contribution to the thermal 
conductivity extracted using (||) to calculate the to- 
tal quasiparticle contribution Kf^{B,0). Note that the 
field dependence of tana|| is affected by the value of 
obtained from (^, in particular the linear dependence 
with the magnetic field observed above '--^ 2 T. Decreas- 
ing 5% the obtained kP|J the computed tanan decreases 
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FIG. 9: (a) Left axis: Magnetic field dependence of tan ay ob- 
tained using the phonon thermal conductivity wg^ obtained 
from the best fit of k,xx{B) using (^ (•), or decreasing this 
by 5% (o). Right axis: Magnetic field dependence of the ra- 
tio /3. The dashed and straight lines are guides to the eye. 
(b) The same as panel (a) but the points were obtained from 
Eq. (^ assuming pure d^2_y2-wave symmetry and with the 
same parameters that fit the field dependence of the oscilla- 
tion amplitudes in Fig. ^. 



substantially at high fields, see Fig. H This result em- 
phasizes the difficulty to obtain true field dependences of 
thermal transport parameters by using phenomenological 
separation methods. 

A magnitude free of the separation model is given 
by the ratio between the averaged change with field of 
the longitudinal thermal conductivity, {kxx{B = 0) — 
i^xx{B, 6))g and the averaged transverse thermal conduc- 
tivity {\Kxy{B)\) e ^ which only contains information of the 
quasiparticle contributions. We define then. 



{MB)\)e 

<Uo,o)-<Uo,s) K'ixiB) 



"xy 



2<(S) 



(10) 



As shown in Fig. ^ (right axis), the ratio (3 increases 
roughly proportional with field. Note first that the sim- 
ple subtraction {n^xiB = 0) — Kxx{B,9))o cancels the 
phononic contribution of the thermal coHdjuetiijity which 
we assume it does not change with fieldd'ElliSEJ. Second, 
(3 does not provide the same information as the ratio 
Kxx/'^xy Therefore, /3 can be used for comparing with 
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any model for the quasiparticle transport. We note that 
the observed Unear regime of (3 should disappear with 
increasing field. When the field is increased, the vortices 
are moving closer together and hence the scattering prob- 
ability of QP by vortices increases. Thus, at higher fields 
than those used in the present work, either a saturation 
or a non-monotonic behavior in the field dependence of 
as well as is expected. 
With Eq. (||) and with the same parameters that fit 
the oscillation amplitudes (see Fig. 0) we have calculated 
both parameters tan(a||) and /3, see Fig. |(b). We found 
that the field dependence of both parameters is qualita- 
tively reproduced by the AS mechanism. However, both 
their absolute values and their slopes with field are not 
given by Eq. (|^) correctly. A disagreement is also found 
if we want to fit the curves given in Fig. fusing the same 
parameters as in Fig. 0. Clearly, the actual AS model is 
not enough to understand the experimental data quanti- 
tatively. 

IV. SUMMARY 

In summary, we have measured the components of the 
thermal conductivity tensor at 13.8 K in an optimally 
doped twinned single crystal of YBCO as a function of 
the magnetic field and the angle 9 between magnetic field 
and heat current. The field-angle dependence measure- 
ments show the predominance of the (1^2 _y2-w&ve pairing 
symmetry over a complex component idxy , which is esti- 
mated to be less than ~ 35% of the resulting total gap, 
within the AS picture for the interaction of QP with vor- 
tices. 

The magnetic field dependence of the amplitude coeffi- 
cients of the fourfold and twofold oscillations observed in 
Kxx and Kxy respectively, allows us to define magnitudes 
free of contributions from the phonons and, hence, suit- 
able for comparison with theory. In particular, the de- 
pendence of the amplitude coefficients n'^y and k^^. , their 
ratio, a Hall-like angle tan a\\ and the ratio (3 between the 
decrease of the longitudinal thermal conductivity and the 
transverse thermal conductivity with field were obtained. 
Assuming only the AS mechanism along with an impurity 
scattering rate, we get a good quantitative description of 
both field dependences of the oscillation amplitude of the 



conductivities only if we increase the intervortex spacing 
five times (or increase the-parallel critical field to 9500 T) 
the value defined in Ref.a. This implies a much smaller 
scattering rate due to AS than the one obtained from the 
original theory. Also, the AS model does not fit with the 
same parameters all the experimental data, although it 
provides a good qualitative description of the field and 
angle dependences of all the parameters. 

The less developed DS model does not allow a rig- 
orous quantitative comparison. It is worth to mention, 
however, that the DS model alone also provides field de- 
pendences of the oscillation amplitudes similar to those 
obtained experimentally. An impurity scattering rate 
should be introduced in the DS theory within a more 
realistic limit than the superclcan limit to perform a de- 
tailed discussion of the quantities involved in the equa- 
tions. The overall results suggest that both AS and DS 
effects are responsible for the observed magnetic field de- 
pendencies but with-different regimes of predominance. 
As suggested in Ref.EZi it is more likely that the DS plays 
an important role at low fields providing a slower mag- 
netic field dependence in the components of the thermal 
conductivity tensor than the AS. This can be understood 
taking into account that while the QP always suffer a DS 
in the mixed state, the probability for Andreev reflec- 
tion increases with increasing field. This might explain 
the less abrupt decrease with field of the ratio n^y/n^xx 
below ~ 2 T and the observed quantitative discrepan- 
cies between the Hall-like angle tana||, the parameter (3 
and the AS model. Further development of the theory 
involving both mechanisms is necessary to confirm this 
hypothesis. 
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